Solar flows and their effect on frequencies of acoustic modes 

Piyali Chatterjee^ 
piyalicStif r . res . in 
and 
H. M. Antia 
antiaStif r . res . in 

Department of Astronomy and Astrophysics, Tata Institute of Fundamental Research, 

Homi Bhabha Road, Mumbai 400005, India 

Received ; accepted 



- 2 - 



ABSTRACT 



We have calculated the effects of large scale solar flows like the meridional 
circulation, giant convection cells and solar rotation on the helioseismic splitting 
coefficients using quasi-degenerate perturbation theory (QDPT). Our investiga- 
tion reveals that the effect of poloidal flows hke the large scale meridional circu- 
lation are difficult to detect in observational data of the global acoustic modes 
since the frequency shifts are much less than the errors. However, signatures of 
large scale convective flows may be detected if their amplitude is sufficiently large 
by looking for frequency shifts due to nearly degenerate modes coupled by con- 
vection. In this comprehensive study, we attempt to put limits on the magnitude 
of flow velocities in giant cells by comparing the splitting coefficients obtained 
from the QDPT treatment with observational data. 



1. Introduction 

Solar convection is believed to be organized in a variety of spatial and temporal scales 
ranging from granules (size ~ 1 Mm, 0.2 hr lifetime), mesogranules (size ~ 10 Mm, 3 hr 
lifetime), supergranules (size ~ 30 Mm, 1 day lifetime) to the giant cells (size > 100 Mm, 1 
month lifetime). The power spectra of convective velocities show distinct peaks representing 
granules and supergranules but no distinct features at wavenumbers representative of 
mesogranules or giant cells (Wang 1989; Chou et al. 1991; Straus, Deubner, & Fleck 1992; 
Straus & Bonaccini 1997; Hathaway et al. 2000). Numerical simulations of solar convection 
routinely show the existence of mesogranules and giant cells (Miesch et al. 2000; Miesch et 
al. 2008). In this study we will only concentrate on the theoretical frequency shifts in global 
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surface gravity and acoustic modes, denoted f- and p-modes respectively due to interaction 
with large scale flows namely the giant cells, the meridional circulation and rotation. The 
meridional circulation, believed to play an important role in magnetic flux transport by 
Solar dynamo modelers (Choudhuri, Schiissler & Dikpati 1995; Dikpati & Charbonneau 
1999; Chatterjee, Nandy & Choudhuri 2004) was observed at the solar surface by Doppler 
measurements of photospheric lines (Duvall 1979; LaBonte & Howard 1982; Komm, 
Howard & Harvey 1993; Hathaway 1996) and later by using techniques of ring diagram and 
time-distance hehoseismology (Giles et al. 1997; Basu, Antia & Tripathy 1999; Gonzalez 
Hernandez et al. 1999). These studies can only measure the flow velocity in the near surface 
layers. This flow having a maximum surface velocity of 30 m s~^ is apparently poleward all 
the way from the equator to the poles at the Solar surface. Conservation of mass requires 
existence of a return flow advecting mass from the poles to the equator somewhere below 
the surface. The dynamo models assume that the return flow occurs near the base of the 
convection zone, but there is no direct evidence about where the return flow is located. It 
is believed that the only hope of detecting such a return flow is by analyzing the effect of 
the meridional circulation on the global helioseismic acoustic modes (Roth & Stix 2008). 
Their investigation using quasi degenerate perturbation theory showed that theoretical 
p-mode frequency shifts due to meridional circulation would be ~ 0.1 //Hz. Such large 
frequency shifts may not be expected from the observed magnitude of meridional flows. For 
the sake of comparison, the equatorial rotational velocity at the Solar surface is ~ 2000 
m s~^ and gives rise to frequency shifts ~ 450 x m nHz as a linear function of azimuthal 
order m. Further, since the flrst order contribution from meridional flow calculated using 
the degenerate perturbation theory vanishes, one would expect the effect to be even smaller 
than what is suggested by the velocity. Thus it is necessary to reexamine these calculations. 

Quasi degenerate perturbation theory (QDPT) was applied to calculate shifts due 
to flows and asphericity in Solar acoustic frequencies by Lavely & RitzwoUer (1992). We 
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have used the same formulation to calculate the p-mode frequency shifts due to differential 
rotation, meridional circulation and giant cells of convection. In contrast to results of Roth 
& Stix (2008), our results indicate a maximum shift of about 1 nHz for a similar meridional 
circulation encompassing the entire solar convection zone. Traditionally, only degenerate 
perturbation theory (DPT) has been used to calculate the effect of rotation on the p-modes. 
Woodard (1989) and Vorontsov (2007) applied QDPT to calculate mixing of eigenfunctions 
due to rotational coupling of p-modes and concluded that the expected line profiles of the 
composite modes in the Solar power spectra are not distorted significantly. In this work, 
we use a different approach to examine if the use of quasi degenerate perturbation theory 
introduces significant corrections in the frequency shifts over that obtained from degenerate 
perturbation theory. It is necessary to check this as inversion techniques for calculating 
rotation rate in the solar interior are based on degenerate perturbation theory. Helioseismic 
inversion for rotation is also insensitive to the N-S asymmetric component of rotation. This 
is because in the degenerate perturbation theory, this effect also gives zero contribution. 
The N-S antisymmetric component of rotation has been studied near the surface and is 
found to be small. In this work we wish to examine its effect to check if its contribution to 
frequency shift are indeed small. 

The giant convective cells have always been elusive to observations at the solar surface. 
In the past there have been studies which have failed to detect giant cell motions (LaBonte, 
Howard & Oilman 1981; Snodgrass & Howard 1984; Chiang, Petro & Foukal 1987) as well 
as those hinting at their existence (Cram, Durney & Guenther 1983; Hathaway et al. 1996; 
Simon & Strous 1997). With the availability of Dopplergrams from SOHO/MDI, it became 
possible to study such long-lived and large scale features more reliably. Beck, Duvall & 
Scherrer (1998) were able to detect giant cells at the solar surface with large aspect ratio 
(~ 4) and velocities ~ 10 m s'^ using the MDI data. Hathaway et al. (2000) used spherical 
harmonic spectra from full disk measurements to detect long-lived power at / < 64. 
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This means that the spatial extent of the features were ~ 2ttRq/^1{1 + 1) > 70 Mm. 
Later Ulrich (2001) identified these long-hved features not as giant cells but higher order 
components of the torsional oscillation. Interestingly it was earlier proposed by Snodgrass 
& Wilson (1987) that torsional oscillations occur due to modulation of differential rotation 
by the Coriolis force arising from meridional motions of giant cells. 

The elusiveness of the giant cell features in the surface velocities lead modelers to 
speculate about the subsurface nature of the giant cell motions (Latour, Toomre & Zahn 
1981). Roth & Stix (1999, 2003) used QDPT to calculate the effect of giant cells on 
p-modes and claimed that giant cells could be found by modeling the asymmetries and 
line broadening in the Solar power spectrum. They claimed that finite line width of 
the multiplets would limit the detection of the frequency splittings to vertical velocity 
amplitude of 100 m s~^ or larger. However, Roth, Howe & Komm (2002) claimed that giant 
cells could be detected with current inversion methods of global hehoseismology as long 
as they exceed an amplitude of 10 m s"^. Usual inversion procedure for rotation neglects 
giant cells and assumes that the odd splitting coefficients arise only from rotation. If there 
are additional contributions to these coefficients, rotation inversions would not give correct 
results. According to Roth, Howe & Komm (2002) the effect of giant cells would appear 
as distortions in the rotation inversion if their velocities are > 10 m s"^. However, using 
the same technique of QDPT we find that with the present global helioseismic data we can 
put an upper hmit of 50-100 m s~^ on the vertical velocity associated with giant cell flows. 
This is mainly because our calculations show much smaller effects from giant cells. 

This paper is organized as follows. In §2 we deflne zonal and poloidal flows inside 
the Sun, while §3 discusses the quasi degenerate perturbation theory, its differences from 

degenerate theory as well as the details of our calculations. Wc present our theoretical 
results and compare them with present f- and p-mode data from both GONG and 
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SOHO/MDI instruments in §4. Finally conclusions are described in §5. 

2. Zonal and Poloidal Solar Flows 

Following earlier works (Lavely & RitzwoUer 1992; Roth, Howe & Komm 2002) we 
express the velocity field in terms of spherical harmonics. For completeness, we give the 
expression here again. 

w;*(r)f X V,F;(^,0)]. (1) 

The quantities m*, v\ and w\ determine the radial profiles of the fiows and V/i is the 
horizontal gradient operator. The Re refers to using only the real part of the spherical 
harmonics as in, 

r [y-*(^,</>) + n*(e,</))]/2 if t is even, 

MY',{eA)] = { (2) 
l[F-*(0,0)-n*(^,0)]/2 if t is odd. 

The first two terms in equation ([T]) define the poloidal component of the fiow whereas the 
last term is the toroidal component. By the poloidal component, we imply the meridional 
and non-zonal toroidal fiows (average over direction is zero) e.g., (i) the meridional 
circulation which carries mass poleward near the surface and sinks near the poles and (ii) 
the giant convection cells, respectively. These fiows are also called large scale fiows to 
distinguish them from other small scale fiows like the turbulent eddies which are of the size 
smaller than the typical scale of global modes used in helioseismology. In presence of only 
the poloidal fiow {w\ = 0) we can apply the equation of mass conservation V.(pov) = to 
get a relation between ul{r) and f*(r) e.g., 

^ 1 ^ 
r or 



sis + l) 



(3) 
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Here po{r) is the density in a spherically symmetric solar model. So now it only remains to 
choose ul{r) appropriately and vl{r) will be determined by equation We choose the 
radial profile of u*(r) as given by equation (19) of Roth, Howe & Komm (2002): 



Here and R^op define the boundaries of region where the flow is confined. The tangential 
component f * also vanishes outside these boundaries. Further, because of very small density 
scale height near the solar surface, f * increases very rapidly near the top boundary and its 
maximum value depends on the choice of -Rtop- For the meridional flow such a behavior has 
not been seen in the actual observed profile (e.g., Basu, Antia & Tripathy 1999). In most of 
our computations we use = 0.7Rq, the approximate position of the base of the convection 
zone, while i?top is taken to be the top of the solar model used in this work, which is around 
1.OO17-R0. With this choice most of the steep variation in f* occurs above r = Rq and the 
variation is probably more reasonable below the solar surface. 

The radial profile of the vertical velocity and the horizontal velocity for the upper 
boundary at i^^p = Rq is shown in Figure [H We wish to point out that Roth, Howe & 
Komm (2002) have used an upper boundary at -Rtop = 0.99/?© in equation (jlj). This avoids 
the region with small density scale height and hence the maximum value of vl{r) is smaller. 
In this work we have used a standard solar model with the OPAL equation of state (Rogers 
& Nayfonov 2002) and OPAL opacities (Iglesias & Rogers 1996) and use the formulation 
due to Canuto & Mazzitelli (1991) to calculate the convective flux. The flow obtained 
using equation (jll) consists of only one cell in the radial direction. Following Roth, Howe & 
Komm (2002) we also use giant cells with a horizontal wave number s = 8. The streamlines 
of the flow for s = 8, t = 8 as well as for s = 8,t = are given in Figure 7 of Roth, Howe & 
Komm (2002). 
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The differential rotation constitutes tlie zonal toroidal flows with t = in last term 
of equation ([T]). The radial and angular dependence of uniform rotation is given by 
Wi{r)dgY^{6,(f)) and that of differential rotation with s = 3 is given by w^{r)dgY.^{6,(j)) 
respectively. Rotation affects the p-modes in two ways: (i) Coriolis force which is linear 
in rotation rate, Q and (ii) centrifugal force which is quadratic in Q. To first order, the 
linear term gives only odd splitting coefficients, while the centrifugal term gives only even 
splitting coefficients. The effect of centrifugal force is of second order compared to (i) and 
has been calculated by Antia et al. (2000). In this work we would consider only the linear 
term due to Coriolis force, or advection in an inertial frame, though we apply QDPT to 
calculate second order corrections due to rotation. 

3. Perturbation of f- and p-modes 

The eigenfunctions and eigenfrequencies for a standard solar model can be calculated 
using perturbative or variational methods. We follow the approach of Lavely & RitzwoUer 
(1992) and use the quasi-degenerate perturbation theory, as opposed to degenerate 
perturbation theory, to determine these quantities. For meridional flow this is necessary 
as the degenerate perturbation gives zero frequency shifts. Let the eigenfunctions of the 
operator Cq governing the seismic oscillations in a standard solar model be called the 
non-rotating standard solar model (NRSSM) eigenfunctions. We shall see in §3.1 how to 
appropriately choose the set of NRSSM eigenfunctions for quasi-degenerate and degenerate 
treatments. 
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3.1. Quasi-degenerate perturbation theory 



The difference between QDPT and DPT lies in the choice of the set of NRSSM 
eigenfunctions (s^), whose linear combinations are used to express the perturbed 
eigenfunctions {s\J e.g., 



k'eK 

In DPT only the eigenfunctions from the standard solar model which are exactly 
degenerate constitute the eigenspace K. In contrast to this in QDPT the set consists of 
all eigenfunctions of the standard solar model with nearly degenerate frequencies. It is 
well known that simple harmonic oscillators couple strongly only if the natural uncoupled 
frequencies of the oscillators are nearly degenerate. Hence it is important to determine which 
of the NRSSM eigenfunctions should be included in the analysis. Let the eigenspace of the 
suitable eigenfunctions be called K. Qualitatively, the quasi-degenerate condition implies 
that we include an eigenfunction Snim{= ini{r)Y{^r + 7]ni{r)VhyD e K iS \ujnim - ^^-refl < ^' 
where e is the radius of the neighborhood about the central frequency. We have generally 
used e = 100 //Hz in this work. This is comparable to asymptotic spacing between two 
consecutive modes in radial order n. Moreover no significant change in calculated frequency 
shifts is observed for e > 100 //Hz. Here cUy-^^ is the frequency of the central eigenmode 
for which we require the frequency shift. The size of the eigenspace K will be determined 
by the (i) desired level of accuracy which is set by e and (ii) the selection rules of the 
perturbation operator. 

In the NRSSM the frequencies are independent of the azimuthal order m and this 
degeneracy is lifted by departures from spherical symmetry due to rotation or magnetic 
field. Helioseismic data traditionally give the sphtting coefficients for all modes that are 
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detected. These coefficients are defined by (e.g., RitzwoUer & Lavely 1991) 

c^n«™ = cu„i + ^4"')pj(m), (6) 

where Uni is the mean frequency of the muhiplet, Vg{m) are the orthogonal polynomials of 
degree q and the so called splitting coefficients. 

The equations of motion for a mode k with eigenfrequency Uk for a NRSSM and 
a model perturbed by addition of differential rotation and/or large scale flow can be 
respectively represented by, 

>CoSk = -poivlsk, (7) 
^os'k + ^^is'k = -Po^'k^^L (8) 

where cu^ is the perturbed frequency and sj^ given by equation (jS]) is the perturbed 
eigenfunction. Taking scaler product with Sj in equation ([H]) and using the notation 
'Hjk' = — f Sj^CiSk'dV and the definition Cis^ = — 2ico'rcfPo(v.V)sk, we obtain the matrix 
eigenvalue equation, 

l^jfe' + h'jiujl - u^^^) I ak' = {uj[^ - ^Jef)ai, (9) 

k'&K 

with, eigenvalue A = {uj'^ — ^^^g^) and eigenvector Xj = {aj}. Here 00^.^^ is a reference 
frequency which approximates u;^. In this work we use uj^.^^ = ook, the frequency of the 
mode being perturbed. It will be clear from §3.2 that the perturbation matrix [H^'^in] is 
Hermitian for both differential rotation and poloidal flows. For example, if we consider only 
two modes with frequencies Ui and UJ2, then quasi degenerate perturbation theory gives the 
following coupling matrix for the mode 002, 



Hn — A H12 

T^21 Ti.22 



(10) 
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with A = cul — cij^. We are interested in the eigenvalues A of [Hij] corresponding to 
perturbation in uj2 which may be easily shown to be, 

^ ~ ^22 - T7 T7 T • (11) 

/til — rt22 — A 

Let be the modified frequency of the mode which initially has a frequency 0^2 • Then 
from equation (fTTl) we have, 



, ,2 1^121 



\2 



(12) 



Till — ^^22 — A_ 

If DPT were used then the last term involving off-diagonal term 7Yi2 will not be present 
and hence this term gives the correction arising from using the QDPT. It is clear that this 
correction can be large if the two modes are nearly degenerate. 

If \Hi2\ ~ |A|, then the frequency shift due to use of QDPT would be of order of A/a;2, 
which is comparable to the difference UJ2 — uJi- On the other hand if |i^i2| <^ |A|, then the 
frequency shift would be much less. For a typical p-mode the spacing uj2 — uJi is of order 
of a few ^Hz. For a rotation velocity of 2000 m s~^, the frequency shift is of the order 
of 450m nHz, which would define the magnitude of diagonal terms. Thus for meridional 
flow with averaged velocity of 20 m s~^ we may expect the matrix elements H\2li02 ~ 4.5Z 
nHz. However, in this case the diagonal elements are zero and we need the off-diagonal 
element which involves cross product of two different eigenf unctions in the integration (cf., 
Eq. [T5l) and (for the same velocity) this integral would be more than order of magnitude 
smaller than that in diagonal term for rotation. Thus the magnitude of ll\2 would be an 
order of magnitude less, giving a frequency shift of order of 0.2/^/ {uj2 — cji) nHz, taking the 
maximum value of m. For / 100 this would give a shift of order of a few nHz. This is 
much smaller than that obtained by Roth & Stix (2008). 

In the next section we shall give the expression for the perturbation matrix elements 
['^ml^ir] "^^^ rotation as well as poloidal flows to be used for solving the eigenvalue 
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equation ([9]) 



3.2. Calculating the perturbation matrix 

The Wigner-Eckart theorem (equation 5.4.1 of Edmonds 1960) states that the general 
matrix element of any tensor perturbation operator can be expanded in terms of Wigner 3j 
symbols whose coefficients of expansion are independent of azimuthal order m and m'. 

^II^S/ = (-ir' I ^ ' ^ \{n'mi\\nl). (13) 
— m' t m 

The coefficient of the Wigner 3j symbol in equation f|T3l) due to coupling by the 
differential rotation is given by, 



(n7'||£*||n/) = Muo^^^{l - (-1)'+^'+^) 



X 



lilv I jsW^spr'^drTs{r), (14) 





where, 



r 



i^'^[/(/ + i) + r(r + i)-.(. + i)]|x 



r^r, I I , (15) 

-1 1 



Tq = \/ q{q + l)/2 and 7g = a/ (2g + l)/47r. ,^ and r/ are radial and horizontal components 
of the eigenfunction defined in §3.1. 

Similarly we have derived the coefficient of the Wigner 3j symbol in equation (fT3!) due 
to coupling by the large scale poloidal flow, 

(n'r||£*||n/) = 87rzcJrgf(l + (-1)'+''+^) x 
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where, 



Re 

lili' I -fsUsPr^dr 
lo 



dr 



r s I 







fdrj dr]' 
dr 



V 



+ 



I' s I 



-1 1 



(16) 



(17) 



Hs{r) 



+ l)-/'(/' + l)] X 
/' s I 


/' s I 
-1 1 



/' s I 
-1 1 



1 -1 



(18) 



The selection rules for the Wigner 3j symbol to be non zero are, 

1- p/2 > max(Z, s), where p = / + + s, 

2. m - m' + i = 0. 

The selection rule for non trivial reduced matrix elements due to differential rotation 
is given by, 

3. l + l' + s must be odd. 



This means that the deY^ component of differential rotation couples modes with I' — I, I ±2. 
The set of modes that constitutes the eigenspace K must follow selection rules (1), (2) and 
(3) and the quasi-degenerate condition \ujnim ~ '^ref I ^ //Hz. The selection rule for 
reduced matrix elements due to large scale poloidal flow to be non zero is 
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4. I + I' + s must be even. 

In this case the set of modes in the eigenspace K must follow selection rules (1), (2) and (4) 
and the quasi-degenerate condition as before. Note that the matrix elements for rotation 
are real and symmetric whereas those for poloidal flow are imaginary and antisymmetric. 
The perturbation matrix \Hn'n'i'i\ hermitian for both cases and will have real eigenvalues. 

4. Results 

4.1. Differential rotation 

We show in this section that QD treatment would give non-zero even coefficients a2q 
even in absence of the centrifugal term and their effect on a2q+i is indeed small as argued 
by Lavely & Ritzwoller (1992). The determination of even splitting coefficients due to 
differential rotation is important as this may be used to remove the effect of rotation from 
the observed splitting coefficients to isolate the effect of magnetic field and other large 
scale flows in the observed splitting coefficients. However, if DPT is used, only the odd 
coefficients a2g+i are non-zero if we neglect the centrifugal force which is second order in Vt. 
Further, in this approximation only the North-South symmetric component of rotation gives 
non-zero frequency shifts. Lavely & Ritzwoller (1992) applied the QDPT to differential 
rotation to flnd that for intermediate I modes (Z ~ 50) the quasi degenerate coupling will 
have little effect on the modal frequencies and can be ignored for all practical purposes. 

It is easy to verify from equation (fT3l) that for pure rotation (v = f2(r, ^)r sin 6'0) the 
symmetric matrix elements are odd functions of m. Hence from equation (I6l) and (|TT1) one 
can infer that the perturbed frequency will have both odd and even splitting coefficients. It 
can be easily shown that there will be no additional effect in QDPT due to component w° 
of rotation velocity as the off-diagonal elements in the resulting matrix would vanish. Thus 
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the leading effect of rotation arises from term. We calculate this contribution by using 



This value is chosen such that the DPT gives a value of close to the observed value for 
most modes trapped in the convection zone. The QDP treatment modifies the odd splitting 
coefficient lai by up to 0.4 nHz (see Fig. [2^) and la^ by up to 0.2 nHz, which are a couple 
of orders of magnitude smaller than the typical errors in these coefficients. In Figure |2)d 
we plot the splitting coefficient la2 for all the modes with frequency less than 4.5 mHz and 
/ > 10, against r"', the lower turning point of the mode. 

The calculation has been done using an equation similar to equation flTTl) but for all the 
modes in the neighborhood of radius 100 /iHz about the central mode. The mathematical 
details of the calculation has been discussed in §3.2 (see equations (fl^ & (fT5l) ). The 
maximum value of la2 ~ 30 nHz which may be significant depending upon the magnitude 
of the large scale flow or magnetic field perturbations. Note that the effect of centrifugal 
force on la2 also happens to be of the same order (see Fig. 5 of Antia, Chitre & Thompson 
2000). We have emphasized earlier that in order to spot the signatures of the large scale 
flows and magnetic fields we need to remove the effect of rotation on a2- 

It has been found from Doppler measurements that there exists a hemispheric 
asymmetry in surface rotation having an angular dependence dgY2 (e.g., Hathaway et 
al. 1996). This component has an amplitude of = —7.8 ± 0.3 m s^^, which indicates 
that the southern hemisphere was rotating slightly faster than the northern hemisphere 
during the period of the data. We have also calculated the splitting coefficient la2 due to 
north-south asymmetry assuming it is constant throughout the convection zone and they 
happens to be small with la2 < 0.1 nHz. 




34.9 m s"^ if r > 0.7R 



otherwise. 



(19) 
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4.2. 



Meridional circulation 



In this section we investigate theoretical frequency shifts due to meridional circulation. 
We start with s = 2,t = 0, which is the dominant component in the observed meridional 
flow near the surface. In our analysis, we have included all f- and p-modes with frequency 
less than 4.5 mHz and / > 10. We consider a meridional circulation with one cell in each 
hemisphere and maximum horizontal velocity at the surface equal to 30 m s^^ which is 
consistent with the Doppler measurements and ring diagram analysis. This gives Uq = 9 m 
in equation (jlj). From equations f|T71) & (|T8l) we see that the coupling matrix for poloidal 
flows is Hermitian with zero diagonal elements. This means that degenerate perturbation 
treatment cannot to be applied to calculate the frequency shifts and it becomes essential to 
use the QDPT, unlike for rotation where QDPT just provides a second order correction. In 
presence of zero diagonal elements of the matrix [Hij] , equation (ITT!) for the frequency shift 
for coupling between two modes simplifles to, 



depends on the sign of A = uj'^^^ — ujf, where ui is the nearest mode to the central frequency. 
If there are many modes with frequency close to frequency of the central mode, then we can 
expect these contributions to be added. In this case, if the sign of frequency differences are 
not the same, the terms will partially cancel each other. If for some pair of modes which 
satisfy the selection rules, the frequency difference is very small, then their frequencies 
would be shifted signiflcantly. 

In Figure [3] we show an example of the frequency shift 5z/(m) of the multiplet 
(n, /) = (1,292), which was also considered by Roth & Stix (2008). The Wigner 3j symbol 
with t = (implying m' = m) in equation f|T3l) is an odd function of m and so 5z/(m) oc Ti^g 
is symmetric about m = for all the multiplets. Another important point to note from 




47rco'2A 



(20) 



It may be noted that the sign of the frequency shift for the central multiplet with uJ-^^f 
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equation (l20l) is that 5v oc Uq. This has also been verified by numerical calculations. The 
maximum shift we get for this multiplet from our calculations is ~ 8.5 nHz in contrast to 
Fig. 2 of Roth & Stix (2008) where they obtain a maximum shift of 0.1 yuHz. On comparison 
with their equation (15) we find a difference of a factor of 2 in the first term in the expression 
of Hs{r). Even after using their version of equation (15) in our calculations we were not 
able to reproduce the large shifts reported by them. The reason for this discrepancy is not 
clear. One possibility is normalization of velocity amplitude. The velocity decreases rapidly 
with height and if the normalization is applied at higher level, then the entire velocity 
profile will be scaled up. In these calculations we have normalized mq by comparing it with 
the surface velocity. We find that the peak value of f*(r) in Figure [1] occurs a little below 
the surface. For instance for a fiow with s = 2 and a velocity of 30 m s~^ near r = Rq, 
this peak in v\{r) happens to be 72 m s~^ at r = O.996i?0. However, in reality such a rapid 
increase in Vs hasn't been found from local helioseismology which can probe up to a depth 
of O.96i?0. Hence, in our opinion the theoretical frequency shift calculated using such a 
radial profile of v\{r) is an upper limit to the actual shift. Ideally, we should normalize uq 
to get the maximum value of the horizontal velocity to be 30 m s~^. That will bring down 
the splitting coefficients by a factor of 4 or more. On the other hand, if we normalize the 
velocity further up then for the same value of velocity the calculated splitting coefficients 
would go up. However, that is not realistic as in that case the maximum velocity would be 
much larger than the observed value. In Figure |lK,b we plot the frequency shifts averaged 
over m for each multiplet where as Figure shows the splitting coefficient la2 calculated 
for all the multiplets. 

We consider only even values of s so that the meridional fiow does not have cross 
equatorial components. As we increase s, selection rules allow more and more multiplets in 
a radius of 100 yuHz to couple and there are many more instances of 'near degeneracy', so 
that the 5z/ oc 1/A in equation fl20l) for some multiplets becomes quite large. For example 



- 18 - 



the multiplets (17, 56) and (16, 64) which have a frequency difference of —0.06 /iHz can be 
coupled by the meridional flow with s = 8 to give la2 ~ 480 nHz. In these calculations we 
have used Uq = 100 m s^^. In Figure [3] we have plotted the frequency shifts for the pair of 
multiplets (15, 34) and (14,40) also considered by Roth, Howe & Komm (2002), coupled by 
the same flow. The result for all multiplets with frequency less than 4.5 mHz and / > 10 is 
shown in Figure [S^,b. It is to be noted from Figure \5jp that only some multiplets have 02 
larger than the errors in 02 for those modes in observational data. It is only these multiplets 
which we may hope to detect in the observations. In this case we have taken the errors from 
the GONG data set centered at 19 November 2002. We present comparisons with GONG as 
well as MDI data in §4.4. Apart from the nearly degenerate modes the splitting coefficients 
of other modes are larger than that found for s = 2 with similar value of horizontal velocity 
near the surface. It may be noted that for large values of s the maximum value of horizontal 
velocity would be quite different from f*(r) because of the additional factor arising from 
the gradient of spherical harmonics (cf., Eq. [1]). In the case where we normalize uo such 
that v^IRq) = 30 m s~^, we also obtain an increase with s in the splitting coefficients as 
reported by Roth & Stix (2008). This increase is not seen if the maximum value of the 
horizontal velocity near the surface is normalized to 30 m s~^. Of course, observations near 
the solar surface do not show such large magnitudes for these higher order components of 
meridional flow. If realistic amplitudes are used then the effect will be negligible. 



4.3. Giant convection cells 

In this section we calculate the splittings due to poloidal flows with angular dependences 
that depend on longitude, 0. We consider the flow given by spherical harmonics Y^{6, 0) 
and Yg{6,(j)) (also called sectoral rolls or banana rolls). This calculation is little more 
involved than that in §4.2 as a non zero t allows coupling between different m and m' in 
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equation (|T3l) . Hence it becomes important to take into account the effect of rotational 
splitting on A before calculating the effect of these kind of poloidal flows. It is easy to see 
from the properties of the Wigner 3j symbols that the s = 8,t = 8 flow couples the mode 
(n, /, m) with (n', /', m ± 8). Hence the difference of the square of frequencies A is no longer 
independent of m unlike in §4.2. 

Let us consider the perturbation in the frequency of the multiplet (15, 34). The nearest 
multiplet to this happens to be (14,40) with a frequency difference of 4.12 /iHz and the 
next one is (16, 28) with frequency difference 19.01 yuHz. This is the same mode shown 
by Figure 3 of Roth, Howe & Komm (2002). We use the rotational splittings calculated 
from temporally averaged rotation rate as inferred from the GONG data to calculate 
A. The result of our calculation including only the nearest multiplet (14,40) has been 
shown in Figure El The asymmetric frequency shift shown by the solid line in Figure ^ 
denotes the change in the frequency of the mode (15, 34, m) due to coupling with the 
mode (14, 40,m — 8) whereas the dashed line denotes the change due to coupling with 
(14, 40,m + 8). As shown in equation ( l20i) a crucial component of the shift comes due 
to the difference of the squared frequencies A^,™' of the modes coupling according to 
the QDPT condition and the selection rules. The corresponding Am,m-8 and Am,m+8 are 
plotted in Figure Eb- The total shift is given by the sum of the solid and dashed curves. 
The couplings within a multiplet i.e., {n, I, m) ^ (n, /, m ± 8) happen to be zero because of 
the anti-symmetry of the matrix elements. The values of splitting coefficients lai and la2 
for the mode (15, 34) are calculated to be —12.2 nHz and —5.7 nHz and that for (14, 40) are 
11.1 nHz and 5.3 nHz, respectively. 

We repeat this calculation for all modes with frequency less than 4.5 mHz and find 
that some of the multiplets have quite high values of 02 (see Fig. (Hb). For example, 
the multiplet (18, 61) which couples with nearest multiplet (17, 69) with a frequency 
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difference of —3.4 /iHz. In Figure [7^ we have plotted 5v{m) versus m for the couphng 
(18, 61, m) (17, 69, m — 8) which seems to have a discontinuity at m = —22. This is due to 
the fact that „j_8 for this couphng also changes sign at m = —22. Note that two modes 
(18,61,m) and (17, 69,m — 8) having a mean frequency difference of 4.12 /iHz become 
'nearly' degenerate at m = —22. This jump presumably gives rise to high values of splitting 
coefficients. The other coupling (18, 61, m) ^ (17, 69, m + 8) is smooth and so is the A^.m+s 
as shown in Figure [Tjo by the dashed line. We also calculate the total shift experienced by 
the multiplet (18,61) while coupling with 3 nearest modes (17, 69), (19, 55), (19, 53) and 
find no significant difference with Figure [Th-. We have also shown a calculation for the 
asymmetric shifts for the multiplet (14,34) and (14,40) coupled due to a flow varying as 
Yi{e, 0) in Figure El 

The asymmetric shift due to these couplings give rise to non-zero odd splitting 
coefficients in addition to the even coefficients. Since rotation inversions neglect this 
contribution, the result may not be correct if the giant cells have significant amplitudes. 
The effect of giant cells on rotation inversions can be estimated by performing rotation 
inversion using odd splitting coefficients from these calculations. The resulting rotation 
rate should be subtracted from the actual rotation inversions. It is likely that if these 
multiplets are used for helioseismic inversion using the rotation kernel, they will give rise 
to distortion in the inverted rotation profile. The effect on the odd splitting coefficients, 
Oi and as is shown in Figure [9K,c in terms of observational errors in those multiplets. We 
use the theoretically calculated ai and to perform a 1.5d helioseismic rotation inversion 
using the formulation due to Ritzwoller & Lavely (1991). We use 1.5d inversion as we have 
calculated only 2 splitting coefficients, ai and 03. For this purpose, we use only those modes 
which are present in a GONG data set and use the errors in those modes for inversions 
using a Regularized Least Squares method using the same smoothing as used for inverting 
real data (e.g., Antia et al. 1998). The results are shown in Figure [TOl which shows the 
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resulting Wi{r) and w^lr) in terms of the estimated errors. Because of presence of modes 
with large 01,03 with ~ 0.7Rq, there are some oscillations in this region, but their 
amplitudes are less than the estimated errors. Thus the effect on rotation inversion may 
not be significant even when the giant cells have velocity of 100 m s~^. With increasing 
velocity the effect would be larger and it may be possible to detect this effect as it will give 
an oscillatory signal in rotation inversion at these depths. This has been pointed out by 
Roth, Howe & Komm (2002). They find a much larger effect as their estimate of frequency 
shift is generally larger than what we find. In real data the effect may not be as dramatic 
as shown by Figure 8 of Roth, Howe & Komm (2002) since many of the modes with large 
shifts in 02^+1 would yield large residuals in inversions, and any reasonable inversion would 
eliminate such modes. In our calculation, we found residuals going to about 4(T level and 
did not eliminate any modes. Thus it is difficult to put any limits on magnitude of flows 
using inversion results. 

The even coefficient 02/^2 is also shown in Figure [9)3. If we find some multiplets with 
0.6 < ti/Rq < 0.8 having magnitudes larger than the errors in observations as in Figure [9b 
we should be able to throw light on the flow velocities in giant cells. We shall look for such 
features in 54.4. 



4.4. Comparison with observed data 

We use 7 different data sets from the GONG observations, covering the late descending 
part of cycle 22 to end of cycle 23, i.e., from June 1995 to November 2008. We choose 
only those modes for which 02 is available in the data sets. Since the different data sets 
can have different set of modes we repeat the calculation for all the 7 data sets. Most of 
the contribution to observed splitting coefficients appear to come from near surface effects, 
while we are interested in modes with turning point close to the base of the convection 
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zone. Thus we first separate the surface effects in the data by fitting a cubic spline in terms 
of frequency to llni0'2k/ ^nioQik versus frequency. Here, Ini is the mode inertia of the mode 
(e.g., Christensen-Dalsgaard 2002), while Inio is the mode inertia for the / = mode at the 
same frequency, while Qik is a geometric factor as defined by Antia et al. (2001). After 
subtraction of the surface effect from the observed splitting coefficients, all the multiplets 
are sorted in the order of increasing and then an error weighted 100 point average is 
applied. In Figure [TTk . b we show the result of the averaging for the GONG 108-day 
averaged data sets centered on 19 November 2002 and 18 September 2007 respectively. To 
get an idea of the magnitude of the giant cell flows one can consider comparing the size 
of the valley and the hump in the region 0.6 < ti/Rq < 0.8 in the theoretical curve with 
that in the observation. It is difficult to make out if the expected signature is present in 
the observed data sets. In most cases no significant hump is seen in this region, though in 
some cases like in Figure [TTh there is some hint of hump which is comparable to errorbars. 
The main problem with this approach is that the running mean tends to average modes 
with positive and negative values of 02, thus reducing the significance of the results. Thus 
it may be better to separate out these modes to look for the signal. 

We have said earlier that only way to put an upper limit on the flow velocity is to 
look for signatures of nearly degenerate modes in the data. In order to do that we divide 
the multiplets from the theoretical calculations in two groups: one for which la2 > 10 nHz 
and the other with la2 < —10 nHz. Amongst 3700 theoretically evaluated multiplets we 
find only about 250 multiplets satisfying either criteria. Then we further search for these 
multiplets in the GONG data sets and perform an averaging over the available multiplets 
to get an estimate of 02 from theory and error (T2 in 02 from the observations. Only about 
10-12 of the 250 multiplets are usually found to be present in the GONG data sets. We find 
that for all the data sets, \la2\ ~ 22 nHz for theoretically calculated coefficients, whereas 
the error in la2 (in observed set) is ~ 14 nHz for the Yg{6,(j)) flow with uq = 100 m s~^. 
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The corresponding values for the Y^{9, cf)) flow with a i^o = 50 m s~^ are 28 nHz and 13 
nHz respectively. Since the observed data do not show these values (in both > 10 nHz 
and < —10 nHz sets) we can only put an upper limit on velocities in such flows. For this 
purpose we compare the 02 values from theory with observational errors to calculate the 
confidence level of the upper limit. So we can say that up to a 1.5a level we can rule out 
the Yq{6, (J)) dependent flow with uq > 100 m and the existence of the ^8^(^, 0) flow with 
1^0 > 50 m can be eliminated with a confidence level of 2a. We have also performed 
some calculations with Yg{9, 0) angular dependence and can say with a 1.5a confidence 
level that such fiows also cannot have > 70 m s~^. Our results for multiplets in all 
the 7 GONG data sets divided into two groups according to theoretical value of la2 are 
represented in Table 1. MDI data also give similar results. By looking at detailed peak 
profiles of these abnormal modes in observed data it may be possible to put more stringent 
limits on velocity of such flows. 



5. Conclusions 

In this study we have calculated the effect of differential rotation, the meridional 
circulation and the giant cell fiows on p-mode frequencies using quasi-degenerate 
perturbation theory. For toroidal fiows like differential rotation, the quasi degenerate theory 
provides only a second order correction to the eigenfrequencics, whereas for the poloidal 
flows it becomes necessary to use the quasi degenerate treatment since the diagonal terms 
of the perturbation matrix vanishes. We agree with Lavely & RitzwoUer (1992) that the 
effect of rotation on the odd coefficients is negligible and hence using the degenerate theory 
is sufficient. However, the even splitting coefficients are about one-half of the errors in 
observations and they might not be neghgible while calculating the a2q for other effects hke 
the magnetic field. This contribution is comparable to the effect of centrifugal term. 
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The frequency shift due to N-S antisymmetric component of rotation rate vanishes in 
the degenerate perturbation treatment and it is necessary to use QDPT. We find that with 
reahstic magnitude of velocities in the antisymmetric component, the computed sphtting 
coefficients are very small and this effect may be neglected. 

For a meridional circulation with maximum horizontal velocity at the surface of 30 
m s~^ and one cell per hemisphere we find that la2 ~ 8.5 nHz for multiplets with turning 
points rt near the surface. Because of our choice of velocity field, the horizontal velocity 
increases rapidly with depth near the surface. This behavior is not seen in the Sun. In 
that case it may be more meaningful to choose velocity profile with a maximum value of 
30 m s~^, which occurs a little below the surface. In that case the splitting coefficients 
would need to be reduced by a factor of 4 or more, making them even smaller. The mean 
frequency shift due to meridional flow is found to be up to 12 nHz, which is much smaller 
than what is found by Roth & Stix (2008). The reason for this discrepancy is not clear. 
This could be due to differences in where the velocity is normalized to specified value. 

In any case, la2 is very small for multiplets with rt/ Rq < 0.9 coupled with single celled 
meridional circulation. This value is less than roughly one-seventh of the errors in the 
observed splitting coefficients. We also do not find any change with increasing the number 
of cells in the radial direction or by changing the depth of penetration of the meridional 
circulation while keeping the surface amplitude constant. This makes it impossible to 
comment on the return flow believed to be present at the base of the convection zone. When 
we increase the number of cells per hemisphere the selection rules allow more multiplets to 
couple with each other. In the process some of the nearly degenerate multiplets combine to 
give rather large values of la2. Interestingly, most of these nearly degenerate modes have 
turning points in the range O.Gi?© to O.Si?©. These purely meridional flows give frequency 
shifts which are symmetric about m = and only give non zero values for even splitting 
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coefRcients. The observed amplitude of these higher order components near the surface 
is rather small and with such amplitudes the effect is expected to be small. It appears 
that in general the magnitude of splitting coefficients due to = 9 m s~^ with s = 2 is 
smaller than that with x^o = 100 m s~^ and s = 8, even after excluding modes with nearly 
degenerate frequencies. This could be due to a larger maximum horizontal velocity as well 
as larger radial velocity (for s = 8) or some geometric factors arising in the calculations. 

In addition to this, we have also calculated the asymmetric frequency shifts caused 
due to giant convection cells with angular flow profiles given by Y^{6,(j)) and Yg{6,(j)). 
Our results differ significantly from earlier works of Roth, Howe & Komm (2002) and 
Roth & Stix (2008) who find the effect of poloidal flows on p-mode frequencies to be an 
order of magnitude larger than what we flnd. The asymmetric shifts also contribute to the 
theoretically calculated odd coefficients. We have performed a 1.5d rotation inversion on 
these coefficients to detect any discernible feature in the inverted profile. However, since 
the magnitude of the features in the inverted profile are smaller than the inversion errors 
we conclude that giant cells with uq < 100 m s~^ do not have much effect on the rotation 
inversion. 

Finally we have analyzed some of the GONG data sets covering solar cycle 23 to look 
for possible evidence of giant cells. We do not find any signal of these fiows in selected 
modes with large splitting coefficients and we can put upper limits on velocities of such 
cellular fiows. Prom our analysis one can say that the existence of giant convection cells 
of Yg{9,4>) angular pattern with maximum vertical velocities of 100 m s~^ can be ruled 
out with a confidence level of 1.5a whereas cells with an angular dependence ¥^{9, (f)) 
cannot have vertical velocities > 50 m with a confidence level of 2a. It is important to 
remember here that the GONG data sets are averaged over 108 days, while MDI data sets 
over 72 days. This may be somewhat larger than the lifetimes of giant cells and hence the 
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signal may be averaged out. We can look at shorter data sets, but in that case the errors 
would be larger. Further, the convective cells are not expected to have smooth velocity 
profiles of the form used by us over the entire solar surface. That would also reduce the 
contribution to splitting coefficients. 

Finally, we would like to point out that while the degenerate perturbation theory can 
be easily adopted for inversion because of linearity of effect and the form of perturbation 
which can be represented by appropriate kernels. In contrast, the effect of quasi-degenerate 
perturbation theory is nonlinear in velocity and further is mainly governed by frequency 
differences between nearly degenerate modes, rather than the matrix elements which depend 
on velocity profile. Hence it is difficult to use these shifts in inversions. For example, the 
main effect of flow with s = 8 is felt in modes with turning points close to 0.7Rq, which 
happens to be close to the lower limit of region in which our flows were localized. But 
we have verified that this is just a coincidence by performing calculations with different 
lower boundary. The same set of modes is affected irrespective of lower boundary, because 
it is these modes that have small frequency differences. This is only determined by the 
value of s which determines the range of / values that are coupled. Thus various types of 
perturbations with the same coupling will give large splittings in these modes and it is 
difficult to identify the actual source of perturbation. 

We thank the referee for useful comments and suggestions which have helped in 
improving the manuscript. This work utilizes data obtained by the Global Oscillation 
Network Group (GONG) project, managed by the National Solar Observatory, which 
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Foundation. The data were acquired by instruments operated by the Big Bear Solar 
Observatory, High Altitude Observatory, Learmonth Solar Observatory, Udaipur Solar 
Observatory, Instituto de Astrofisico de Canarias, and Cerro Tololo Inter-American 
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Fig. 1. — Radial profiles of vertical velocity Ug{r) (thick solid line) and horizontal velocity 
v^{r) (thin solid line) with uq — 100 m and U2{r) (thick dashed line) and horizontal 
velocity V2{r) (thin dashed line) with uo — 9 m s~^. 
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Fig. 2. — (a) Odd splitting coefficients lai as a function of rt for the dgY^ component of 
differential rotation (refer to equation ([1])). (b) Same as (a) but for even splitting coefficient 
la2. 
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Fig. 3. — Frequency shifts as a function of azimuthal order m/l. The shift in the multiplet 
{n, I) — (1, 292) due to couphng by Y2 flow with = 9 m (sohd hne); (n, V) — (15, 34) 
due to flow and = 100 m (dashed line); shift for the multiplet (n, X) — (14, 40) for 
the same flow and velocity (thick solid hne) . 
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Fig. 4. — Mean frequency shifts (5i/(m) as a function of (a) frequency u and (b) turning point 
Tt/ Rq due to coupling by meridional flow with s — 2,t — (c) Splitting coefficient la2 as a 
function of lower turning point radius r^. 



-35- 



-20 









* * 












* 






4 0.5 0.6 0.7 0.8 0.9 1 




• 

• 


• 










• 



• 










• 
• 

• 






4 0.5 0.6 0.7 0.8 0.9 1 

r / R 
t s 



-1 



Fig. 5. — (a) Mean frequency shift Si/ as a function of turning point rt/ Rq due to coupling 
by large scale flow with s = 8, i = 0. (b) Splitting coefficient 02 as a function of turning 
point radius r*. The points denote the theoretical values of 02/0-2' for modes (n,/). The 



errors have been chosen from the GONG data set centered around 25 December 2002. 
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Fig. 6. — (a) The solid line gives the frequency shift because of interaction between 
modes (15, 34, m) and (14, 40, m — 8) due to the flow with an angular dependence Y^{9, 4>). 
The dashed line is the frequency shift due to interaction between modes (15,34,m) and 
(14,40,m + 8). The total shift in frequency for the mode (15,34) is the sum of the solid 
and the dashed lines. The frequency shift for the mode (14, 40) is given by the dotted line. 
(b) A^^^_8/2i/ (solid line) for the coupling (15,34, m) ^ (14, 40, m — 8); and Am,m+8/'^i^ 
(dashed line) for the couphng (15, 34, m) ^ (14, 40, m + S). These values are in /xHz. 
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Fig. 7. — (a) The solid line gives the frequency shift because of interaction between 
modes (18, 61, m) and (17, 69, m — 8) due to the flow with an angular dependence Yg(9, (f)). 
The dashed line is the frequency shift due to interaction between modes (18,61,m) and 
(17, 69,m + 8). The total shift in frequency for the mode (18,61) is the sum of the solid 
and the dashed hues. The frequency shift for the mode (17, 69) is given by the dotted line. 
(b) A„^m_8/2z/ {solid line) for the coupling (18, 61,m) ^ (17, 69, m — 8); and Am,m+8/2i^ 
{dashed line) for the coupling (18, 61, m) ^ (17, 69, m + 8). These values are in /xHz. 
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Fig. 8. — Asymmetric shift in frequency of the multiplet (a) (15, 34) (sohd hne) upon couphng 
with neighboring multiplets (14,40) and (16,28) in presence of large scale flow s — — 
4. The dashed-dotted line is the shift for the mode (14,40) due to coupling with (15,34) 
and (13,48). (b) Same as (a) but for (18,61) (solid line) upon coupling with neighboring 
multiplets (17, 69) and (19, 55) in presence of large scale flow s = 8, i = 4. The dashed-dotted 
line is the corresponding shift in (17, 69) due to coupling with 2 nearest neighbors. 
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Fig. 9. — (a) tti/ai and (b) 02/(72 and (c) 03/(73 as a function of lower turning point radius 
rt for the Yg(9,(f)) kind of flow with uq — 100 m s~^. The (71,2,3 are the errors in the 
corresponding observational splitting coefficients from the GONG data set centered about 
November 2002 (see Table. 1). 
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Fig. 10. — (solid line) and (dotted line) in terms of the respective errors are shown as 
a function of radial distance obtained from the 1.5d inversion of coefficients oi and 03 shown 
in Figure 9. 
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Fig. 11. — la2 as a function of turning point radius after averaging over 100 nearest points 
in Tt for GONG data set centered around (a) November 2002 and (b) 18 September 2007. 
The dotted hne correspond to the Y^{9, 0) flow; dashed to Y^{9, 0); thick sohd to Y^{9, (p). 
The sohd hne with errorbars show the observed data points. To avoid congestion only a few 
representative errorbars are shown. 
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Table 1: Comparison of observed la2 with theory for flows with an angular dependence 
Y^{9,(f)). The numbers in the GONG data set have the following correspondence with the 
date about which the 108-day GONG month is centered: (1) 29 June 1995, (2) 18 June 1997, 
(3) 18 September 2000, (4) 19 November 2002, (5) 23 July 2004, (6) 18 September 2007 and 
(7) 2 July 2008. For i = 4, = 100 m s"^ whereas for i = 8, -Uq = 50 m s'^. Note that for 
each t the group I of multiplets correspond to la^^ < — 10 nHz whereas the group II is for 
la^ > 10 nHz. CL is the confidence level on the upper limit of uo- 



Data t 
set 


1. laf < -10 nHz 11. laf > 10 nHz 


1(^2 laT CL 


7 nViQ 7 7 til /^T 


1 4 
8 


17.7 16.9 -22.6 1.3 
-11.3 14.5 -29.0 1.9 


-3.7 14.8 21.9 1.5 
3.3 13.1 24.2 1.8 


2 4 
8 


-8.7 14.9 -22.2 1.5 
-23.4 13.2 -28.1 2.1 


-12.4 14.0 21.3 1.5 
4.7 11.7 24.8 2.1 


3 4 
8 


-24.5 14.6 -21.7 1.5 
-30.2 13.8 -29.7 2.1 


-21.8 14.3 22.6 1.6 
-17.7 12.3 24.2 2.0 


4 4 

8 


-1.3 14.5 -22.3 1.5 

-30.2 13.8 -29.7 2.1 


-12.0 14.3 21.8 1.5 

-25.2 12.0 23.9 1.9 


5 4 
8 


-28.3 14.2 -22.0 1.5 
5.3 13.1 -27.3 2.0 


-5.8 13.9 22.0 1.5 
-4.4 11.8 24.1 2.0 


6 4 
8 


-0.4 12.3 -22.3 1.8 
-0.6 10.7 -28.1 2.6 


-22.1 12.8 22.2 1.7 
-4.5 10.4 26.0 2.5 


7 4 
8 


-13.2 14.2 -22.2 1.6 
-27.4 12.7 -28.9 2.2 


-27.0 13.5 21.8 1.6 
-0.6 11.4 24.6 2.2 



